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ABSTRACT 

We obtained thermal equilibrium solutions for optically thin, two-temperature black hole 
accretion disks incorporating magnetic fields. The main objective of this study is to explain 
the bright/hard state observed during the bright/slow transition of galactic black hole candi- 
dates. We assume that the energy transfer from ions to electrons occurs via Coulomb collisions. 
Bremsstrahlung, synchrotron, and inverse Compton scattering are considered as the radiative 
cooling processes. In order to complete the set of basic equations, we specify the magnetic flux 
advection rate instead of j3 = J? g as/Pmag- We find magnetically supported (low-/?), thermally 
stable solutions. In these solutions, the total amount of the heating via the dissipation of turbu- 
lent magnetic fields goes into electrons and balances the radiative cooling. The low-/? solutions 
extend to high mass accretion rates (> a 2 A/Edd) and the electron temperature is moderately 
cool (T e ~ 10 s — 10 9 5 K). High luminosities (> O.liEdd) and moderately high energy cutoffs in 
the X-ray spectrum (~ 50 — 200 keV) observed in the bright/hard state can be explained by the 
low-/? solutions. 

Subject headings: accretion, accretion disks — black hole physics — magnetic field — X-rays: binaries 



1. Introduction 

Galactic black hole candidates (BHCs) arc 
known to exhibit transitions between different X- 
ray spectral states. Typically, a transient outburst 
begins in the low/hard state at a low luminos- 
ity. The X-ray spectrum in the low/hard state 
is roughly described by a hard power law with a 
high energy cutoff at ~ 200 keV. As the luminos- 
ity increases, these systems undergo a transition to 
the high/soft state (so-called a hard-to-soft tran- 
sition). The X-ray spectrum in the high/soft state 
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is dominated by the disk emission of characteristic 
temperature ~ 1 keV. 

Recently, two distinct types of hard-to-soft 
transitions, the bright/slow transiti on and the 



dark/fast transition, are reported (e.g.. lBelloni et al.l 
l2006tlGierlinski fc Newtonhood ) . The bright /slow 
transition occurs at ~ 0.3 LEdd and takes more 
than 30 days. The system undergoes a transition 
from the low/hard state to the high/soft state via 
the "bright/hard" state and the very high/steep 
power law (VH/SPL) state during the bright/slow 
transition. The X-ray spectrum in the bright/hard 
state is described by a hard power-law with a 
(moderately) high energy cutoff at ~ 50 — 200 keV 
and the luminosity is "brighter" than that in the 
low/hard state. The dark/fast transition occurs at 
less than 0.1 LEdd and takes less than 15 days. The 
system immediately switches from the low/hard 
state to the high/soft state during the dark/fast 
transition. Here, LEdd = iircGM/ k cs ~ 1.47 x 
10 39 (M/10M©) (k os /0.34 cm 2 g-^erg s" 1 is 
the Eddington luminosity, M is the black hole 
mass, and k cs is the electron scattering opacity. 
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Mivakawa et al. ( 20081 ) analyzed the results of 
RXTE observations of the BHC GX 339-4 in 
the rising phases of the transient outbursts (this 
object showed the bright/slow transition in the 
2002/2003 outburst and the dark/fast transition 
in the 2004 outburst). They found that the cutoff 
energy strongly anti-correlates with the luminos- 
ity and decreases from ~ 200 keV to ~ 50 keV 
in the bright /hard state, while the cutoff energy 
is roughly constant at ~ 200 keV in the low/hard 
state. This suggests that the electron temperature 
of an accretion disk emitting hard X-rays decreases 
as the luminosity increases in the bright/hard 
state. Furthermore, the bright/hard state has 
been observed from ~ 0.07 ^Edd up to ~ 0.3 Lsdd- 
They concluded that such anti-correlation is ex- 
plained by the scenario that the heating rate from 
protons to electrons via the Coulomb collision bal- 
ances the radiative cooling rate of inverse Comp- 
ton scattering. The main purpose of this paper 
is to present a model explaining the bright/hard 
state. 

In the conventional theory of accretion disks, 
the concept of phenomenological a- viscosity is in- 
troduced. In this framework, the -nu^-component 
of the stress tensor, which appears in the an- 
gular momentum equation and the viscous heat- 
ing term, is assumed to be proportional to the 
gas pressure, — — assPgas (we ignore the ra- 
diation pressure in this paper because we focus 
on optically thin disks). He re ass is the viscos- 
i ty pa rameter introduced by IShakura fc Sunvaev 
(|l973l ). The magnetic field can be an origin of 
the a-viscosity because the Maxwell stress gener- 
ated by the magneto-rotational instability (MRI) 
efficiently transp orts angular momentum in accre- 



tion disks (e.g., iBalbus fc Hawlevl 119911 ) and the 



dissipation of magnetic energy co ntributes the disk 
heating (e.g.. iHirose et al.l l2006). 

Optically thin, hot accretion disks have been 
studied to explain h ard X-rays from BHCs. 



Thorne fe Pried (|1975l) proposed that hard X-rays 
from Cyg X-I a re produced in an inner op tically 



thin hot disk. Shibazaki & Hoshi ( 1975) stud 



opt 

m 



ied the structure and stability of optically thin 
hot accretion disks. We have to consider two- 
temperature plasma in such disks because the 
electron temperature is expected to become lower 
than the ion temperature in such a low density, 
high temperature region. The energy equations 



for ions and electrons are written in the form 
Pi T i ^ = (l-6 heat )q + -q ie , 



PeT e 



dS 
~~dt 



5heatg + + q lc - q Iad 



(1) 



(2) 



where q + is the viscous heating rate, g le is the 
energy transfer rate from ions to electrons via 
Coulomb collisions, q^ is the radiative cooling 
rate, and <5heat is the fraction of heating to elec- 
trons. The left hand-sides represent the heat ad- 
vection terms for ions (<?ad,i) and electrons (<7ad,e)- 
We note that early works on optically thin, hot, 
two-temperature accretion disks assumed that the 
viscous heating acts primarily on ions (5hcat <1). 
Eardlev et all (jl975l ) and IShapiro et ail (|l976l 



hereafter SLE) constructed a model for optically 
thin two-temperature accretion disks. In the SLE 
solutions, the dissipated energy is transferred from 
ions to electrons via Coulomb collisions (q + ~ q ic ) 
and radiated away by electrons (q le ~ q^d)- Al- 
though the electron temperature is high enough to 
explain the X-ray spectrum in the low/hard state, 
the SLE solutions are thermally unstable in the 
framework of the a-prescription of viscosity. 



Ichimarul ([19771 ) pointed out the importance of 



heat advection in hot, magnetized accretion flows, 
and obtained steady solutions of optically thin 
disks. Such geometrically thick, optically thin, 
advection-dominated accretion flows (ADAFs) 
or radiatively inefficient accretion flows (RIAFs) 
have b een s tudie d extensively by Naravan fc Yil 
(|1994L Il995l) and lAbramowicz et all f|l995fT In 
the ADAF/RIAF solutions, a substantial frac- 
tion of the dissipated energy is stored in the 
gas as entropy and advected into the central 
object (q + ~ <7ad,0- Only a small fraction of 
the dissipated energy is transferred to electrons 
and radiated away. Th e ADAF/R I AF so lutions 



are thermally stable. Esin et all (1997) found 



that the maximum mass accretion rate of the 
ADAF/RIAF solutions is M c ,a ~ l-3a 2 Msdd 
which corresponds to L ~ 0.4a 2 LEddj where Msdd 



i s the Eddington mass accretion rate. lEsin et al 



(Il998h showed that the electron temperature in 
the ADAF/RIAF solutions weakly anti-correlates 
with the luminosity and decreases to ~ 10 9 5 K. 
These features are consistent with the facts that 
the energy cutoff weakly anti-correlates with the 
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luminosity around ~ 200 keV in the low/hard 
state, and that these systems undergo a transi- 
tion from the low/hard state to other X-ray spec- 
tral states (i.e., the high/soft state during the 
dark/fast transition, and the bright/hard state 
during the bright/slow transition) about at this 
maximum luminosity of the ADAF/RIAF solu- 
tions. However, the ADAF/RIAF solutions can- 
not explain the strong anti-correlation in the range 
of L > 0.1 L E dd and T c < 200 keV observed in the 
bright/hard state. 

The heat advection works as an effective coolin 
in the ADAF/RIAF solutions. lYuanl (|200ll [200 
presented a luminous hot accretion flow (LHAF) 
in which the heat advection for ions works as an ef- 
fective heating. Above the maximum mass accre- 
tion rate of the ADAF/RIAF solutions, the heat 
advection overwhelms the viscous heating and bal- 
ances the energy transfer from ions to electrons 
(<7ad,i ~ q lc )- The LHAF s are t hermally unsta- 
ble. However, I Yuan et"aL ( 2003 ) concluded that 
the thermal instability will have no effect on the 
dynamics of the LHAFs because the accretion 
timescale is shorter than the timescale of growth of 
the local perturbation at such high mass accretion 
rate. The LHAF solutions also cannot explain the 
bright/hard state because the electron tempera- 
ture is high and roughly constant at ~ 10 9 ' 5 K. 

In these models mentioned above, magnetic 
fields are not considered explicitly, and the ra- 
tio of the gas pressure to the magnetic pres- 
sure is assumed to be constan t (typi cally, /? — 
Pgas/Pmag ~ 1)- IShibata et al. (1990) suggested 
that an accretion disk evolves toward two types 
of disks, a high-/3 disk and a low-/? disk, by 
carrying out two-dimensional magnetohydrody- 
namic (MHD) simulations of the buoyant escape 
of the magnet i c flux owing to the Parker insta- 
bility (jParkerl Il966h . In the high- ,5 disk, the 
magnetic flux escapes from the disk owing to the 
Parker instability and /? inside the disk is main- 
tained at a high value. Global three-dimensional 
MHD simulations of optically thin, radiatively 
inefficient accretion disks also indicate that the 
amplification of magnetic fields becomes satu- 
rated when /? ~ 10 in a quasi-steady state ex- 
cept in the pl unging region very close to the 
black hole (e.g..lHawlevll200filMachida et alfc oOO; 



once a disk is dominated by the magnetic pres- 
sure, it can stay in the low-/? state because the 
strong magnetic tension suppresses the growth of 
the Parker instability. 

iMachida et al.l ([20061) demonstrated that an op- 
tically thin, radiatively inefficient, hot, high-/? 
(ADAF/RIAF- like) disk undergoes transition to 
an optically thin, radiatively efficient, cool, low-/? 
disk except in the plunging region when the mass 
accretion rate exceeds the threshold for the onset 
of a cooling instability. During this transition, the 
magnetic flux (B V )H is almost conserved at each 
radius because the cooling timescale is shorter 
than that of the buoyant escape of the magnetic 
flux, where (B v ) is the mean azimuthal magnetic 
field and H is the half thickness of the disk. In this 
way, the magnetic pressure becomes dominant and 
supports the disk as the gas pressure decreases ow- 
ing to the cooling instability. Eventually, the disk 
stays in a quasi-steady, cool, low-/? state. Because 
the MRI is not yet stabilized in this quasi-steady 
state, the magnetic field still remains turbulent 
and dominated by the azimuthal component. As 
a result, the heating owing to the dissipation of 
the turbulent magnetic field balances the radia- 
tive cooling. 

Johansen fc Levin performed local three- 



Hawlcv fc Krolikl 2001 



20031 IMachida et al.l 120041 ) . On the other hand. 



Machida fc Matsumotol 



dimensional MHD simulations of strongly mag- 
netized, vertically stratified accretion disks in a 
Keplerian potential. They showed that strongly 
magnetized state is maintained near the equatorial 
plane because the buoyantly escaping magnetic 
flux is replenished by stretching of a radial field. 
The MRI feeds off both vertical and azimuthal 
fields and drives turbulence. The Maxwell and 
Reynolds stresses generated by the turbulence be- 
come significant. Therefore, they indicated that 
highly magnetized disks are astrophysically viable. 

We note that such low-/? disks are quite differ- 
ent from magnetically dominated accretion flows 
(MDAFs; iMeierl [2005) observed in the plunging 
region of optically thin accretion disks in global 
MHD and general relativis tic MHD simulations 
(e.g., Fragile fc Meier 20091 ) in terms of their en- 
ergy balance and configuration of magnetic fields. 
Outside the plunging region, the magnetic field be- 
comes turbulent and dominated by the azimuthal 
component because the growth timescale of the 
MRI is shorter than the inflow timescale. As a re- 
sult, the released gravitational energy is efficiently 
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converted into the thermal energy via the dissi- 
pation of the turbulent magnetic field. On the 
other hand in the plunging region, the ratio of the 
timescales is reversed because the inflow velocity 
increases with decreasing the radius in such disks. 
Therefore, magnetic field lines are stretched out 
in the radial direction before turbulence is gener- 
ated by the MRI and is dissipated. As a result, 
a substantial fraction of the gravitational energy 
is converted into the radial infall kinetic energy. 
Since the heating owing to the dissipation of tur- 
bulent magnetic fields becomes inefficient, the gas 
pressure becomes low, and the flow becomes mag- 
netically dominated. Although both of the low-/3 
disk and the MDAF are cool and magnetic pres- 
sure dominant, they are essentially different. We 
focus on the low-/3 disk in this paper. 

Oda et ail ( 2007 ) constructed a steady model 



of optically thin, one-temperature accretion disks 
incorporating magnetic fields on the basis of these 
results of three- dimensiona l MHD simulations of 
accretion disks. Oda et all ( 2009() extended it to 
the optically thick regime. They assumed that the 
t£7</3-component of the stress tensor is proportional 
to the total pressure. In order to complete the 
set of basic equations, they specified the advec- 
tion rate of the azimuthal magnetic flux instead 
of (3. They found a new thermally stable solu- 
tion, a l ow-/3 solution, which can explain the re- 
sults by iMachida et~atl <|2006h . In the low- ,3 so- 
lutions, the magnetic heating enhanced by the 
strong magnetic pressure balances the radiative 
cooling. The disk temperature is lower than that 
in the ADAF/RIAF solutions and strongly anti- 
correlates with the mass accretion rate. They also 
found that the low-/? solutions exist above the 
maximum mass accretion rate of the ADAF/RIAF 
solutions. Therefore, they concluded that the opti- 
cally thin low-/3 disk can qualitatively explain the 
bright/hard state. However, they considered one- 
temperature plasma and bremsstrahlung emission 
as a radiative cooling mechanism. It is expected 
that the electron temperature becomes lower than 
the ion temperature and that synchrotron cooling 
and/or inverse Compton scattering become effec- 
tive in such disks. 

In this paper, we extend the model of opti- 
cally thin, one-temperature disk to that of opti- 
cally thin, two-temperature disks. We consider 
synchrotron emission and inverse Compton scat- 



tering as a cooling mechanism in addition to 
bremsstrahlung emission. We obtained the ther- 
mal equilibrium curves and found that the opti- 
cally thin low-/3 solutions can quantitatively ex- 
plain the bright/hard state. The basic equations 
are presented in Section[2] In Section[3l we present 
the thermal equilibrium curves. Section |4] is de- 
voted to a discussion. We summarize the paper in 
Section El 

2. Models and Assumptions 
2.1. Basic Equations 

We extended the basic equations for one- 
dimensional steady, optically thi n, two-temperatu re 



black hole accretion flows (e.g.. iKato et al.l 12008) 
incorporating magnetic fields. We adopt cylin- 
drical coordinates (zu, ip, z). General relativistic 
effects are simulated using th e pseudo-Newtonian 
poten tial ip — —GM/(r — r s ) (jPaczvnskv fc Wiita 
Il980l ). where G is the gravitational constant, M 
is the black hole mass (we assume M = 10M Q in 
this paper), r = (vj 2 + z 2 ) 1 / 2 , and r s = 2GM/c 2 is 
the Schwarzschild radius. For simplicity, the gas 
is assumed to consist of protons (ions) and elec- 
trons. The number density of ions and electrons 
are equal by charge neutrality, n = ni = n e . 
We start with the resistive MHD equations 



dp 

at 



9v , _. 
W + („.V)„ 



V • (pv) = , 



-pVip - Vpg 



(3) 



jxB 



, (4) 



dt 



d (p c e c ) 

at 



+ v • [(m+pi)v] - (v-v)pi 



(l-S hcat )q+ -q ic 



+ V • [(p e e e + p e ) v]-(v- V) p e 



<5heat9 + + q m - q. 



rad 



(5) 



(0) 



— =VX^ XJ B-- W j , (7) 

where p = p; + p c is the density, pi = m\n and 
p c = m c n are the ion and electron densities, m-^ 
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and m c are the ion and electron masses, v is the 
velocity, B is the magnetic field, j = cV x B/An 
is the current density, p gas = Pi+p c = nk (2] + T c ) 
is the gas pressure, pi and p c are the ion and elec- 
tron gas pressure, T\ and T e are the ion and elec- 
tron temperature, k is the Boltzmann constant, 

£ i = (Pi/pi) / (7i ~ 1) and e c = (p e /A>) / (7e - 1) 
are the internal energy of ions and electrons. Here, 
7i = 5/3 and j c — 7 C (T c ) are the specific heat 
ratio for ions and electrons. In the energy equa- 
tions for ions ([5]) and electrons ([6]), q + is the heat- 
ing rate, g~ d is the radiative cooling rate, and q m 
is the energy transfer rate from ions to electrons 
via Coulomb collisions. Here, <5heat represents the 
fraction of heating to electrons. In the induction 
equation ([7]), rj m = c 2 /An a c is the magnetic diffu- 
sivity, where a c is the electric conductivity. 



dv v dvu, 

PVvj-^ h PVz-z— 

dm dz 



m 



J__d_ 

m 2 dm 



2 \ B vj B v) 
An 



d_ 

dz 



[B V B Z 
An 



_ _chp_ _ 1 dptot 
dz p dz 



(10) 
(11) 



d v d 

^— + Pi) v w ] + — (piQ + Pi) + — [(piei + pi) v z ] 
aw w oz 

9 9 ft x ^+ 

dm oz 



(12) 



d v d 

— [(p e e e + p e ) v m ] + — {p c e c + p c ) + — [(p c e c + p e ) v z ] 
OVJ VJ oz 



9 9 A 4 

dm dz 



2.1.1. Azimuthally Averaged Equations 

Three-dimensional global MHD and local radiation- 
MHD simulations of black hole accretion disks 
showed that magnetic fields inside the disk are 
turbulent and dominated by the azimuthal com- 



ponent in a quasi-steady state fe.g. jMachida et al 



20061 Iffirose et alll2006f ). On the basis of results 
of the simulations, we decomposed the magnetic 
fields into the mean fields B = (0,B V ,6) and 
fluctuating fields SB = (SB^^SB^^SBz) and also 
decomposed the velocity into the mean veloc- 



ity V : 
ity Sv 



, v z ) and the fluctuating veloc- 



(Sv w , Sv v , Sv z ). We assume that the 



fluctuating components vanish when azimuthally 
averaged, (Sv) = {SB) = 0, and that the radial 
and vertical components of the magnetic fields are 
negligible compared with that of the azimuthal 
component, \B V + SB V \ ^> \SB Z \. Here ( ) 

denotes the azimuthal average. 

Let us derive the azimuthally averaged equa- 
tions. We assume that the disk is in a steady 
state and in hydrostatic balance in the vertical 
direction. By azimuthally averaging equations ([3]) 
- ((7|) and ignoring the second order terms of Sv, 
SB m , and 5B Z , we obtain 



1 d 

7{ (mpv z 
m dm 



^-(pv z ) = 
dz 



(8) 



PVv 



' dm 



+pv z 



dz 



PK 



= -P 



dil> d Ptot (B%) 
dm dm 



Anm ' 
(9) 



+{V x (Sv x SB)} V - {rj m V x (V x §)} v , (14) 

where ptot = Pgas + Pmag is the total pressure and 
Pmag = (Bp)/8ir is the azimuthally averaged mag- 
netic pressure. The third and fourth terms on the 
right-hand side of Equation (fT4"]l represent the dy- 
namo term and the magnetic diffusion term which 
we approximate later on the basis of the results of 
the numerical simulations. 

2.1.2. Vertically Integrated, Azimuthally Aver- 
aged Equations 

We assume that the radial velocity v w , the 
specific angular momentum I — mv v , and /3 = 
Pgas/Pmag are independent of z, and that the disks 
are isothermal in the vertical direction. Under 
these assumptions, the surface density E, the ver- 
tically integrated total pressure W to t , and the half 
thickness of the disk H are defined as 



pdz 



Pa exp 



2H 2 



dz = V2np H , (15) 



Ptotdz 



/'...inCXp j ~^jpj ) '/"- = \ 2-/) : „,aMij 



(17) 
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where f2 K o = (GM/m) 1/2 /(m - r s ) is the Keple- 
rian angular velocity. Here the subscript refers 
to quantities in the equatorial plane. Using the 
equation of state for the ideal gas, the vertically 
integrated total pressure is expressed as 



W, 



gas 



mag 



^(l + r 1 ) .(is) 



Now we integrate the other basic equations in 
the vertical direction. We obtain 



(19) 



dv^ 1 dWtot 

V^n— h 



dm £ dm 
£ 2 - l|p W tot d\nn K0 2/3- 1 Wtot 1 



dm 



1 + /3- 1 £ 



,(20) 



M(£- 4 



-27TH7 2 ^ {B ^ BtP/ 



Q 



ad . ; 



(i 



-'heat 



47T 



, (21) 



(22) 



/ad 



,e = <WtQ + + Q 1C - Q 



rad ' 



(23) 



v m {By)dz 



[{V x (Sv x JB)}^ 
- {ry m V x (V x B) } v ]dmdz + const(24) 

where M is the mass accretion rate, £ko — m 2 flKo 
is the Keplerian angular momentum and £[ n is 
the specific angular momentum swallowed by the 
black hole. In the energy equations, Q ac i,i and 
Q a d,e are the vertically integrated heat advection 
terms for ions and electrons, Q + , Q~ ad , and Q lc 
are the vertically integrated heating rate, radia- 
tive cooling rate, and energy transfer rate from 
ions to electrons via Coulomb collisions. In Equa- 
tion <I> is the radial advection rate of the 
azimuthal magnetic flux (hereafter we call it the 
magnetic flux advection rate). 

In this paper, we assume that £ = £ko instead of 
Equation (|2"D|) because we focus on local thermal 
equilibrium solutions. 



2.2. ^-Prescription of the Maxwell Stress 
Tensor 

Global MHD simulatio ns of radiatively ineffi 



cient, accretion flows (e.g., Hawlev fc Krolikll2001 



Machida et al. l2006h showed that the ratio of 
the azimuthally averaged Maxwell stress to the 
sum of the azimuthally averaged gas pressure 
and magnetic pressure is nearly constant («b = 
-(B TU B ip /4:TT)/(p gas + Pmag ) - 0.05 - 0.1) except 
in the plunging region. Following the simulation 
results, we assume that the azimuthally averaged 
tuy'-component of the Maxwell stress inside a disk 
is proportional to the total (gas and magnetic) 
pressure 



47T 



-Otptot ■ 



(25) 



Integrating in the vertical direction, we obtain 



H 



-in 



-dz = — ctWu 



(26) 



This is one of the key assumptions in this paper. 
When the magnetic pressure is high, the stress 
level can be high even though the gas pressure is 
low. 

We can rewrite this relation in terms of the 
kinematic viscosity, v, as 



v = A v a\J c s0 2 + c A0 2 H 



where 



Q, rflnO 
f^KO d\nm 



(27) 



(28) 



CsO 



= v Pgaso / Po is the sound speed, and cao = 



\/2pmago//3o is the Alfven speed. 

2.3. Prescription of the Magnetic Flux 
Advection Rate 

We complete the set of basic equations by spec- 
ifying the radial distribution of the magnetic flux 
advection rate. If we perform the integration in 
the second term of the induction equation (|24l) . 
we obtain 



-v w B (m)V47rH 

[dynamo and diffusion terms] 



(29) 
const. 
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where 



4 f kT[ + fcT e 



1/2 



E l/2^-l/ 2( g-l/ 2 (3Q) 



is the mean azimuthal magnetic field in the 
equatorial plane. According to the result of 
the global three-dim ensional MHD simulation by 
Machida et al. ( 20061) . the magnetic flux advection 
rate at a radius is roughly unchanged before and 
after the transition from the ADAF/RIAF-like 
disk to the low-/? disk. Hence, we adopt the mag- 
netic flux advection rate as the parameter in order 
to complete the set of the basic equations. The 
magnetic flux advection rate depends on various 
mechanisms such as the escape of magnetic fluxes 
due to the magnetic buoyancy, the regeneration 
of azimuthal magnetic fields by the shear motion, 
the generation of magnetic turbulence through 
the MRI, dissipation of magnetic fields due to the 
magnetic diffusivity, and magnetic reconnection. 
If the sum of the dynamo term and the magnetic 
diffusion term is zero in the whole region, the 
magnetic flux advection rate is constant in the ra- 
dial direction. The globa l three-dimensional M HD 
simulation performed by lMachida et al" ( 20061 ) in- 
dicated that the magnetic advection rate increases 
with decreasing radius, specifically, $ oc w' 1 , in 
the quasi steady state as a result of magnetic dy- 
namo and diffusivity processes. Because it is hard 
to compute the dynamo term and the magnetic 
diffusion term explicitly from the local quantities, 
we parameterize the dependence of $ on w by 
introducing a parameter, C,, as follows. 



*(w;C,M) = *aut(M) 



w 

^out 



-C 



(31) 



where ^out is the magnetic flux advection rate at 
the outer boundary w — vo OVlt - When £ = 0, the 
magnetic flux advection rate is constant in the ra- 
dial direction. When £ > 0, the magnetic flux ad- 
vection rate inc reases with decre asing radius (See 
also Figure 1 in lOda et a.l.ll2009h . Here we deter- 
mine the parameter $ ut by imposing the outer 
boundary condition, Ti, ut = T e , ut = ^Virial = 
[(mi + m c ) c 2 /3fc] (tn out /r s ) _1 and /3 out = 10 at 
tj7 out = 1000r s . This leads that $ out oc M 1 / 2 . In 
the local model presented in this paper, the value 
of £ just means the amount of the magnetic flux 
advection rate at a radius (e.g, $(nj = 5r s ) = 



$ out , 4.9 $ out , 24 $ out for C = 0, 0.3, 0.6, respec- 
tively) . 

Equation pip is the second key assumption in 
this paper. Specifying the magnetic flux advec- 
tion rate enables the magnetic pressure to increase 
when the disk temperature decreases. By contrast, 
if we specified the plasma /3 at each radius instead 
of the magnetic flux advection rate, the decrease 
in the temperature results in a decrease in mag- 
netic pressure. This is inconsistent with the re- 
sults of three-dimen sional MHD simulations (e.g., 
Machida et"aTll2006l ). 

We address the similarity between the concepts 
of the mass accretion rate and the magnetic flux 
advection rate in order to facilitate understand- 
ing of the concept of the magnetic flux advection 
rate. The right-hand side of equation (fT9|) repre- 
sents the mass flux crossing at zu per unit time and 
we have denoted it by M . If there is no mass loss 
or gain (e.g, due to inflows and/or outflows), M is 
constant in the radial direction. Otherwise M is 
a function of m (e.g., M oc uj s , s is a parameter). 
The right-hand side of equation (f29| represents 
the azimuthal magnetic flux crossing at w per unit 
time and we have denoted it by $. If there is no 
loss or gain of azimuthal magnetic fields (e.g, due 
to the dynamo and/or the magnetic diffusion), <1> 
is constant in the radial direction. Otherwise $ is 
a function of w and we prescribed it as equation 
PT|) . That is, when £ = 0, the magnetic flux is 
conserved in the radial direction, and when C > 
(or £ < 0), the magnetic flux increases (or de- 
cr eases) with a dec reasing radius (see Figure 1(b) 
in lOda et a.1.1 12009^ . We note that the azimuthal 
magnetic flux inside a disk can increase when the 
azimuthal flux of opposit e polarity buoyantly es- 
capes from the disk (e.g.. iNishikori et al]|2006l ). 



2.4. Energy Equations 

2.4-1- The Magnetic Heating Rate 

In the conventional theory, the viscous heating 
was expressed as qj~ is = t^ v w (dfl/dw) and as- 
sumed to heat primarily ions (^hcat ~ rn c jm\ ~ 
10 -3 ), where is the tz7</?-component of the to- 
tal stress and f2 is the angular velocity. 

The results of three-dimensional MHD simu- 
lations indicate that the dissipation of the mag- 
netic energy dominates the total dissipative heat- 
ing rate throughout a disk and is expressed as 
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q + ~ (B 777 B in /4:Tr)w ( dQ/dv c ) Ce.g.. iHirose et al 
2006; iMachida et aD 120061 : iKrolik et al.l l2007t l 



Hereafter, we refer to it as the magnetic heating 
rate. Following these simulation results, we em- 
ploy magnetic heating as the heating mechanism 
inside a disk, and set the vertically integrated 
heating rate as follows: 



cr 



— ; w ~, — 

47T dm 



dz 



dQ 

dvo 



(32) 



where we have used equation (1231) . We note that 
if the magnetic pressure is high, the heating rate 
can also be high even when the gas pressure is 
low. The eventual expressi on of the heating term 
is not at all unusual (e.g.. iNakamura et al.lll997 ; 
Manmoto et al.lll997t lYuanll200ll ) except that we 
consider much lower values of 0. We will discuss 
in Section 14.11 a lower l imit of below which the 
MRI is stabilized (e.g.. iPessah fc Psaltis1l2005l ) so 
that this expression is no longer valid. 

We note that the magneti c heating d o es not 
always heat primarily ions. I Yuan et al. (2003) 
suggested that <5hcat ~ 0.1 — 0.5 is required to 
fit the spectr u m of Sgr A* with RIAF models. 
Sharma et al.l (|2007l ) carried out local shearing 
box simulations of the nonlinear evolution of the 
MRI in a collisionless plasma considering pres- 
sure anisotropy and showed that 5hoat ~ (1 + 
Sy/Ti/Tc)^ 1 (note that the definition of <Sh oa t = 
at Iq + is different f rom that of <5h oa t = it lit m 
Sharma et al.l 120071 ) . The exact value we choose 
for this parameter is not so important, in partic- 
ular, for the low-/? solutions. We find that the 
low-/? solutions presented in this paper practically 
unchanged for any value of 5heat 0.1. Therefore, 
we adopt <5 hcat = (1 + "i^jTi/T^ 1 as a fiducial 
value. 

2.4-2. The Energy Transfer Rate from Ions to 
Electrons by Coulomb Collisions 

If the ion temperature is higher than the elec- 
tron temperature, Coulomb collisions transfer en- 
ergy from ions to electrons. The energy transfer 
rate from ions to electrons per unit volume via 
Coulo mb collisions is given bv lStepnev fc Guilbert 
(Il983h 



3 m c 
2 m\ 



n 2 o~TC (In A) 



JzTq k'-L \ 



2(0 e + di) 2 + 1 



2K Q 



(33) 



where ctt is the Thomson scattering cross sec- 
tion and In A is the Coulomb logarithm (roughly 
In A ~ 20). K n are modified Bessel function of 
the second kind of the order n, respectively. The 
dimensionless electron and ion temperatures are 
defined by 



kT c 
m c 2 



kT; 



(34) 



For technical reason, we use the following for- 
mula which uses no special functions, and is ac- 
curate to within a factor of 2 when B\ < 0.2 
(jDermer et al.lll99lh 



i0 3m c 2 i.T^ 27r ) 1/2 
q = n orcm A(fcJ — ki\) — 

2 TO; 



Integrating in the vertical direction, we obtain 
3m P S 2 



Q 1C = -o 



2 mi 2(m; + m c ) 2 y/nH 



(Txcln A 



■X 



K 2 (l/e a )K 2 (l/0i' 



{kT e - kTi 



(27T 



,1/2 



(0e 



a/2 



(0 C + 0i) 



a/2 



(35) 



(36) 



> 0.2) 



< 0.2) 



We note that 6*i < 0.2 in almost all solutions 
presented in this paper. 

2.4-3. Radiative Cooling Rate 

We assume that the radiative cooling occurs 
through electrons and consider bremsstrahlung, 
synchrotron, and Compton cooling by bremsstrahlung 
and synchrotron photons as cooling processes. 
The vertically integrated radiative cooling rate 
is expressed as 



Qrad — *5br + Qsy + Qbr,C + ^sy, 



sy,C 



(37) 



Fo l lowing Naravan fc Yi dl995l) (see also Svenssonl 
19821: IStepnev fc Guilbertl ll983). bremsstrahlung 



cooling rate per unit volume is 

%r = Cei + 9br.ee = n 2 (j T ca { m c c 2 [F ci (6 c ) + F ee (9 e )] , (38) 



S 



where the subscripts ei and ee denote the electron- 
ion and electron-electron bremsstrahlung cooling 
rates, at is fine-structure constant, the function 
F c i(9 c ) and F ee (0 o ) have the approximate form 



90c 
2ir 



ln(2?70e + 0.48) 



= < 



2ft 



1/2 



(0 C > 1) 



4 Hf [1 + 1.78W C 134 ] (d e <l) 



emissivity can be approximated by the blackbody 
emission from th e surface of the disk. Following 
lEsin et ah ( 19961 ). we estimate v c as the frequency 
at which the synchrotron emission from the region 
w ~ m + Anj is equal to the blackbody emission 
(in the Rayleigh- Jeans limit) from the upper and 
lower surfaces of the region. This condition gives 
(39)the equation 

v 2 

(2irvjf\vj)E sy dv = 2(2-Kw^m)2n-^kT c dv . (45) 



-*■ Of 



9ft 

7T 
5 



— -(hx(2r]9 e ) + 1.28) 
(44-3tt 2 )^ /2 x 



>D 



6tt 3 / 



[ (i + i.w c + e 2 c - i.25^o /2 ) {e c < 1 



r]B = exp(— 7e) and 7e ~ 0.5772 is Euler's num- 
ber. Integrating in the vertical direction, we ob- 
tain 



: br 



(mi + m c ) 2 y/nH 
[F ci {8 c ) + F CC (9 C )} . (40) 



The synchrotron emissivity of a relativistic 
Maxwellian distribution of electrons in the o p- 
tically thin limit is given bv lPacholczvkl (|l970h : 



e sy dv 



2e n Q 2nv 
V3c K 2 (1/9 C ) 



l'{x M )dv, (41) 



where 



x M 



2v 



eB 
2nm c c 



(42) 



We obtain v c by solving this equation numerically. 
Integrating over frequency, we obtain 



2 / 2-K-^kT e dv ■ 
Jo c 

4irv?kT r 2e 2 



E S ydV 



3c 2 



V3c{m + m e ) if 2 (i/6» ) a } /6 
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+ ^|(a^ c + 3a^ 2/3 



1/3 



fl 4 



0-2 r (19 1/3 



+ 6a 4 t/ c 1/3 + 6) cxp(-a 4 ^ c 1/3 ) 



where the parameters a\, a 2 , 03, and a 4 are defined 
as 



0.4 



0.5316 



«i 



0,2 



3M 2 '^ a\ /i,a3 a{ /2 



-, a 4 = 1 



(46) 



a 1 / 3 , (47) 



with t he fitting function /' (xm) given bv lMahadevan 
(1996) 



et al 



and T(a,x) is the incomplete gamma function 

POO 

T(a,x) = t^ 1 exp(-t)dt. (48) 

J X 

We adopted the prescription for the Comp- 
ton energy e n hance ment factor r\ described by 



Dermer et alj (1991), which is defined to be the 



T , . . 4.0505 / 0.40 0.5316 



1/6 



1/4 



1/2 



average change in energy of a photon between in- 
jection and escape: 



cx P (-1.8899xm /3 ) ■ (43) 
Integrating in the vertical direction, we obtain 



2e 2 



2irv 



-I'{x u )du. (44) 



We assumed that the emission below a critical fre- 
quency, v c , is completely self-absorbed so that the 



V = 1 



where 



P(A-l) 
(1 - PA) 

hv 



1 - 



x 

3ft~ 



-1-lnP/ In A 



(49) 
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p 

m 



1 - exp(-r c . 
PU-l) 



1 - PA 
-1 - In P/ In A 



, A = 1 + 46» 

7?2 = 3 _ '' 3 ?7i, 



1661" 



(50) 



Here, P is the probability that an escaping photon 
is scattered, A is the mean amplification factor in 
the energy of a scattered photon when the scatter- 
ing electrons have a Maxwellian ve locity distribu- 
tion o f temperature 9 C . Following iNaravan fc Yi 
( 19951 ). the vertically integrated Compton cooling 
by bremsstrahlung and synchrotron photons are 
given respectively by 




{ sy,C 



m - m 



»73+l 



(52) 



; ad,< 



M 



2itm 2 rrii + m c 



6 , 



M 



2irzu 2 



TO; 



-Co 



where 
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91nPi 91nS dkiH 

H 

d In vj d In w d In w ' 



(53) 



(54) 



(55) 



Cc = -a e (T e ) 



d In a e (T e ) 



+ 



<91nT e 
<91nE 



d\nt 



dlnT c 
dlnw 
dkiH 
d In vo ' 



(56) 



heat advection term works as an effective cooling 
(heating). 

The entropy gradient parameter for ions, has 
a positive value in ADAF/RIAF solutions, while 
it has a negative value in LHAF solutions. The 
global three-dimensi onal MHD simulations (e.g., 



Machida et al.l 120061 ) and steady, vertically inte- 



?br,C = 3l 7l<9br 



where x c = hv c /m c c?. 

2.4-4- Heat Advection Term 

The vertically integrated heat advection terms 
for ions and electrons are expressed as 



are dimensionless quantities of the order of unity 
(hereafter we call them the entropy gradient pa- 
rameter), en = 1/(1 - 7i ) = 3/2, a c (T c ) = 1/(1 - 
7o(T )), respectively. The positivity (negativity) 
of the entropy gradient parameter means that the 



grated, one-dimensional transonic solutions (e . g . , 
iNakamura et al.l Il996l Il997t lOda et al.l 120071) of 
optically thin black hole accretion flows indicate 
that £i ~ 1. Following these results, we adopt 
£i = 1 as a fiducial value in this paper. We note 
that a value of £i is not important in SLE and 
low-/? solutions because the heat advection term 
is negligible compared to the other terms. 

According to INakamura et al. I dl997l) . £ c can 

(51) 

have a positive and negative value (—0.5 < £ < 
0.5). However, they assumed that no viscous heat 
goes into electrons (<5heat =0). In the energy equa- 
tion for electrons, when <5hcat ^ 0.1, the dissipated 
magnetic energy term is typically greater than or 
comparable to the heat advection term in the in- 
ner region of the disk. Therefore the exact value 
we choose for this parameter is not so important, 
in particular, for SLE and low-/? solutions in which 
<5hcatQ + , Q 1C , Q~ d > <2ad, c - Wc find that the re- 
sults presented in this paper are practically un- 
changed for any value of £ e between —0.5 and 0.5 
when (5h C at 0.1. We show the results for £ c = 0.5 
in most part of this paper. 

3. Results 

We solved the above basic equations at w = 5r s 
for given parameters M, a, (, 5hcat, 6j an d £e- We 
obtained new thermal equilibrium solutions, low- 
13 solutions, in addition to the ADAF/RIAF (for 
positive £j), SLE, LHAF (for negative solutions 
in the optically thin regime. 

3.1. Low-/? solutions 

Figure Q] shows the sequences of each thermal 
equilibrium solution in the S versus M/M-Eddj 
7i (thin line), T e (thick line), and /? plane. The 
disk parameters we adopted are a = 0.05, = 1, 
Co = 0.5, (Shcat = (1 + 3VP7P)- 1 , C = 0.6 (solid), 
0.3 (dashed), and (dotted), respectively. Here 
M E dd = LEdd/VcC 2 = 4ttGM/ (t] c k cs c) is the Ed- 
dington mass accretion rate, rj c — 0.1 is the energy 
conversion efficiency, and k cs — 0.40 cm 2 g _1 is the 
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electron scattering opacity. We obtain three types 
of solutions, ADAF/RIAF (for S < 1 g cm" 2 
at this radius), SLE (for E ~ 1 g cm -2 ), and 
low- /8 solutions (for E > 1 g cm -2 ). We find 
that the low-/? solutions exist above the maxi- 
mum mass accretion rate of the ADAF/RIAF so- 
lutions, M C ,A ~ 0.003AfEdd- This indicates that 
the disk initially staying in the ADAF/RIAF state 
undergoes transition to the low-/? state when the 
mass accretion rate exceeds M C ,A- Furthermore, 
the electron temperature in the low-/? solutions 
is lower (T e ~ 10 8 - 10 9 5 K) than that in the 
ADAF/RIAF solutions. 

The energy balance for ions and electrons is 
illustrated in figure [5] The upper panel shows 
the ratio of the heat advection to the magnetic 
heating for ions, Q a d,i/(1 - Shc&t)Q + (this quan- 
tity is referred to as the advection factor / . 
e.g., Naravan fc Yilll994lll995tlAbramowicz et al 



1995; Yuan 2 00111 . The lower panel shows the ra- 



tio of the heat advection to the total heating for 
electrons (thin line) and the fraction of the energy 
transfer via Coulomb collisions to the total heat- 
ing (thick line). The electrons receive the total 
amount of the magnetic heating in the low-/? solu- 
tions as well as in the SLE solutions even though 
we introduced the parameter <5heat which repre- 
sents the fraction of heating to electrons. The 
fraction (1 — <5heat) of the magnetic heating goes 
into ions and the fraction ^hcat of the magnetic 
heating goes into electrons. However, almost all 
the magnetic heating going into ions is transferred 
to electrons via Coulomb collisions in the low-/? 
solutions ((1 — <5hcat)Q + ~ Q 10 )- Eventually, the 
total amount of the magnetic heating goes into 
electrons. This means that the parameter <5heat 
does not appear practically in the energy balance 
for electrons. The radiative cooling overwhelms 
the heat advection in the low-/? solutions. There- 
fore, the magnetic heating balances the radiative 
cooling in the low-/? solutions (Q + ~ Q^ad)- 

The energy balance is essentially the same in 
both the SLE solutions and the low-/? solutions 
(Q + ~ Qrad)- The difference is which pressure 
dominates the magnetic heating. The magnetic 
heating, which is proportional to the total pressure 
in our model, is dominated by the gas pressure in 
the SLE solutions and the magnetic pressure in 
the low-/? solutions. 

Next we describe the main cooling mechanism 



in each solution. Figure [3] shows the vertically in- 
tegrated bremsstrahlung (solid), bremsstrahlung- 
Compton (dotted), synchrotron (dashed), and 
synchrotron- Compton (dash-dotted) cooling rate 
for ( = (bottom), 0.3 (middle), and 0.6 (top). 
When £ = 0, the synchrotron-Compton cooling 
is dominant in the low-/? solutions for lower mass 
accretion rate while the bremsstrahlung-Compton 
cooling is dominant for high mass accretion rates 
(M > 10 _3 M E dd, E > 1.6 g cm -2 ) even though 
/? < 1. The synchrotron cooling is relatively in- 
effective for higher mass accretion rates because 
of the lower electron temperature. As £ increases, 
the electron temperature become high because the 
large magnetic flux enhances not only the syn- 
chrotron cooling but also the magnetic heating. 
As a result, the synchrotron and synchrotron- 
Compton cooling become efficient. When £ = 0.6, 
the synchrotron-Compton is dominant in whole 
low-/? solutions. 

Now we show that M cx E, T oc E~ 2 , oc E~ 2 
on the low-/? branch (see also Oda et a l. 2009), 
where T = (Tj + T c )/2 is the mean tempera- 
ture. These relations depend on the dependence 
of "JVut on M (our outer boundary condition leads 
$ cx M 1 / 2 ). Here we introduce the parameter 
s, "jout oc M s , in order to leave this depen- 
dence explicitly. First, we derive the relations 
between M and E. Since W to t ~ W mag oc TE/? -1 
on the low-/? branch, equations (fl7|) . ff2T|) - and 
(HP yield H cx T 1 / 2 /? -1 / 2 and M cx TE/? -1 . 
Using equations (JTSJ), (|29"]l. and (|3l]). we find 

that H CX S -(l-2a)/(7-4 S ) ) p K E 2(l-2 s )/(7-4 S ) T) 

and M cx £i-2(i-2*)/(7-4 S )_ Nextj we derive 
the relation between T and E from the energy 
balance of the low-/? solutions (Q + ~ Q~ d ). 
Here we roughly approximate the radiative cool- 



ing rate by Q 



rad 



cx E 2 T 1 / 2 i7- 1 for simplicity 
(this is the same dependence as non-relativistic 
bremsstrahlung cooling for single temperature 
plasma). Since Q + cx Wtot oc M, we find 
that T cx £-2-6(i-2a)/(7-4a) < Therefore /? cx 

S 2(l-2s)/(7-4 s ) T K S -2-4(l-2a)/(7-4») < We find 

that M oc E, T oc E" 2 , and j3 oc E -2 when 
s = 1/2. 

We investigate the relation between the mass 
accretion rate and the electron temperature in or- 
der to explain the anti-correlation between the lu- 
minosity and the cutoff energy observed in the 
bright/hard state. In the low-/? solutions, we ex- 
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pect that the luminosity is roughly proportional 
to the mass accretion rate since Q~ ad ~ Q + oc M, 
and that the electron temperature roughly repre- 
sents the cutoff energy since the inverse-Compton 
scattering is the dominant radiative cooling mech- 
anism. Therefore, the relation between the elec- 
tron temperature and the mass accretion rate 
is useful for comparison with the observational 
data. Figure 0] shows the relations between the 
electron temperature and the mass accretion rate 
for the same parameters as in Figure [TJ We 
find that the electron temperature is typically 
< p]9-5 j£ in the low-/? solutions while it is > 
10 9 - 5 K in the ADAF/RIAF solutions. Further- 
more, the electron temperature in the low-/3 so- 
lutions strongly anti-correlates with the mass ac- 
cretion rate. Since M oc £i-2(i-2s)/(7-4 S )) and 

T OC £-2-6(l-2 s )/(7-4s) in the low _^ somt ions, we 

find that T oc M" 2 " 2 ^ 1 " 2 ^ oc M~ 2 for s = 1/2. 

According to lPessah fc Psaltisl (|2005f ). the MRI 
is stabilized for va > ^/c s vkq (we plot this criti- 
cal point (filled circle in Figure [4j at which va = 
y/c s VKo)- Therefore, the low-/? solutions may not 
exist under the condition that va > y/c s VKo- We 
discuss this issue in Sectional 

We also show the results for a — 0.2 in Fig- 
ure [5] and [6] The maximum mass accretion rate 
of the ADAF/RIAF solution is around M c ,a ~ 
0.05A/Edd- This maximum mass accretion rate 
is higher than that for a — 0.05. Here we in- 
vestigate the dependence of M C ,A on a. Since 
Wtot ~ Wg as , Qad,i/(1 - S hcat )Q + ~ 0.5 and 
Qadx < 5 hcat Q + + Q 1C at M ~ M C ,A, the energy 
equations yield (i) Q ad ,i ~ 0.5(1 - £> h cat)<2 + , (ii) 
Q ic ~ 0.5(1 - <5 h cat)Q + , and (iii)_4catQ + + Q ic ~ 
[(Wt + 0.5(1 - S hcat )]Q + ~ Q rad . Here we as- 
sume <5hcat to be constant for simplicity. We find 
from Equation (i) that T- ~ 10 n K. If Q~ aA can 
be written in a simple form of E 2 /H x f(T c ) like 
Q^ r , we find from Equations (ii) and (iii) that T e 
at M Ci a is independent of a and M c .a is propor- 
tional to a 2 MEdd exactly. Although Q~ ad has a 
complicated form in our model, our numerical re- 
sults indicate that the T c is roughly independent 
of a (T G ~ 10 9 5 K). Using this value of T instead 
of solving Equation (iii), we find from Equation 
(ii) that roughly M C ,A ~ a 2 MEdd- Th is result is 
consistent with M c .a ~ l-3a 2 MEdd bv lEsin et al 



(119971 ) 



3.2. Effect of (5 hG at 

We investigated the dependence of the results 
on <5heat because <5hcat is a poorly constrained pa- 
rameter. We considered three cases: <5h ca t = 0.5 
as an example of model in which ions and elec- 
trons receive equal amounts of the dissipated mag- 
netic energy, <5heat = 10 -3 (~ m c /mi) as an ex- 
ample of a conventional model in which ions re- 
ceive a substantial amount of the energy, and 
<5heat = 0.2 as an intermediate example. We note 
that <5 h cat = (1 + Z^Ti/Te)- 1 lies between ~ 0.05 
and 0.25 in the solutions presented in this paper. 

Figure [7] shows the thermal equilibrium curves 
for the same parameters as in Figure Q] but for 
£hcat = (1 + Z^/WFc)- 1 (solid), 0.5 (dashed), 
0.2 (dotted), and 10~ 3 (dot-dashed), respectively. 
The energy balance for ions and electrons is illus- 
trated in Figure |U 

The ion temperature decreases as <5h ca t in- 
creases because the magnetic heating for ions 
((1 — (5hcat)<3 + ) which is the only heating source 
for ions decreases. M c ,a also decreases because of 
the decrease in the magnetic heating for ions. The 
electron temperature increases as Jheat increases 
in the ADAF/RIAF solutions. By contrast, in the 
low-/? solutions, the electron temperature is almost 
independent of the value of Jheat because 5h oa t does 
not appear practically in the energy equation for 
electrons (6 hcat Q+ + Q ic ~ Q+ ~ Q" d ). There- 
fore, we find that the exact value of <5heat is not so 
important for the low-/3 solutions. 

3.3. Dependence on the entropy gradient 
parameters, £ and £ c 

Figure [9] shows the thermal equilibrium curves 
for the same parameters as in Figure [T] but for 
£i = —1.0 (dashed). We also plotted the thermal 
equilibrium curves for £ = 1.0 (solid) for compar- 
ison. 

We obtained the LHAF solutions in the high 
mass accretion rate and high surface density re- 
gion. The heat advection works as a heating 
for ions and balances the energy transfer via the 
Coulomb collisions in the LHAF solutions. Above 
M C ,A, the heat advection term overwhelms the 
magnetic heating. Hence, the ion temperature be- 
comes higher than that in the ADAF/RIAF solu- 
tions. 

In the LHAF solutions, the electrons receive 
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a larger amount of heat from ions than that in 
the ADAF/RIAF solutions. Nonetheless, the elec- 
tron temperature becomes lower than that in the 
ADAF/RIAF solutions (yet higher than that in 
the low-/? solutions) because the radiative cooling 
becomes more effective in such higher surface den- 
sity region. 

We find that our results are practically un- 
changed for any value of £ e between —0.5 and 0.5 
because the heat advection term is negligible com- 
pared to the other terms in energy equation for 
electrons except when (5h ea t = 10 -3 . Furthermore, 
even when (5h ca t = 10~ 3 , the other quantities ex- 
cept the electron temperature in the ADAF/RIAF 
solutions does not change practically. We show 
the thermal equilibrium curves plotted on the £ 
- T; (thin) and T e (thick) plane for <5h oa t = 10~ 3 , 
Cc = 0.5 (solid), (dashed), -0.5 (dotted), ( = 0.6 
(left panel) , and (right panel) in Figure [101 

4. Discussion 

4.1. Optically Thin, Magnetically Sup- 
ported, Moderately Cool Disks 

We obtained thermal equilibrium solutions 
for an optically thin, two-temperature accretion 
disk incorporating magnetic fields. We included 
bremsstrahlung emission, synchrotron emission, 
and inverse Compton scattering as the radiative 
cooling mechanisms, and introduced the parame- 
ter (Shcat which represents the fraction of the mag- 
netic heating to electrons. We prescribed the zuip- 
component of the azimuthally averaged Maxwell 
stress is proportional to the total pressure. In 
order to complete the set of basic equations, we 
specified the radial distribution of the magnetic 
flux advection rate by introducing a parameter (. 
We found a branch of low-/3 solutions in addition 
to the usual ADAF/RIAF (for & > 0), SLE, and 
LHAF (for & < 0) solutions. 

Here we remark why we can obtain the low- 
/? solutions. First, we prescribed that the zuip- 
component of the Maxwell stress is proportional 
to the total pressure. Therefore, if the magnetic 
pressure is high, we can obtain the magnetic heat- 
ing rate balancing the radiative cooling rate even 
in high surface density and low temperature re- 
gion. Second, we specified the magnetic flux ad- 
vection rate, $, in order to complete the set of 
the basic equations. In the conventional theory, 



j3 is assumed to be constant (typically, (3^1), 
which means that the magnetic pressure is propor- 
tional to the gas pressure. This implies that p gas 
is just multiplied by a constant ((1 — (3~ 1 )p gas ). 
As a result, we cannot obtain the sequence of the 
low-/3 solutions. On the other hand in our model, 
a decrease in temperature results in an increase 
in magnetic pressure (therefore an increase in the 
magnetic heating) under the conservation of the 
magnetic flux in the vertical direction. This is the 
reason why we can obtain the sequence of the low- 
(3 solutions. 

We note that the exact values of £i, £ e) and <5hcat 
are not so important in the low-/? solutions. The 
magnetic heating enhanced by the high magnetic 
pressure balances the radiative cooling in the low- 
(3 solutions. The heat advection terms including 
£i and £ e are negligible for both ions and electrons. 
Furthermore, <5hcat does not appear in the energy 
balance for electrons practically. 

Let us discuss the lower limit of f3. Since 
/? a E -2 and £ oc M under our outer boundary 
condition (s = 1/2), we find that {3 cx M~ 2 , that 
is, /? decreases as the mass accretion rate increases 
in the low-/? solutions. It has been confirmed 
by gl obal MHD simulations (e.g., iMachida et al. 
2006) that the MRI is not stabilized and hence the 
magnetic heating rate is expressed in the form of 
equation ([32]) , at least, when (3 > 0.1. Loc al MHD 



simulations (e.g., Johansen k Levir]|2008l) also in- 
dicated that the Maxwell and Reynolds stresses 
generated by magnetic turbulence are significant 
and yield an effective a- viscosity (a ~ 0.1) in 
highly magnetized disks (/? ~ 1). However, the 
expression of the magnetic heating employed in 
our paper may no longer be valid for much lower 
values of (3 because strong magnetic fields sup- 
press the growth of the MRI. We implicitly as- 
sumed that the dissipation energy of the turbulent 
magnetic fields generated by the MRI is converted 
into the thermal energy of the disk gas. Therefore, 
such heating mech anism becomes ineffective if the 
MRI is stabilized. IPessah k Psaltisl (|2005l) stud- 
ied the evolution of the MRI in differentially rotat- 
ing, magnetized flows beyond the weak- field limit, 
and showed that the MRI is stabilized for toroidal 
Alfven speeds, va = B^/y/Anpo, exceeding the ge- 
ometric mean of the sound speed, c s , and the ro- 
tational speed, vk = mUxo, (i-e., «a > ^/c^vko, or 
equivalently /? < 2c s /wko)- Ovrr results satisfy this 
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condition when f3 < 0.27 (C = 0.6), 0.09 (C = 0.3), 
0.03 (C = 0) for a = 0.05, and /3 < 0.26 (C = 0.6), 
0.08 (C = 0.3), and 0.02 (C = 0) for a = 0.2, re- 
spectively. These critical points are denoted by 
hlled circles in Figure 2] and [5] (see also Table [T]) . 
When /? falls below this critical value, the low-/? 
solutions may not exist because there is no heating- 
source balancing the radiative cooling. As a result, 
the disk may undergo a transition to other states 
(e.g, an MDAF-like disk or an optically thick disk). 

4.2. Thermal Stability 

Let us discuss the thermal stability in this sub- 
section. A general criterion concerning the ther- 
mal instabili t y of disks can b e expressed as (see 
Pringlej[l97l Kato et al.ll2008l ) 



rad/ 



> 



(57) 



In the low-/? solutions, we ignore the heat advec- 
tion term because Q + ~ Q~ d ^> Q a( j. Once 
again we approximate the radiative cooling rate by 
Q~ d cx Y?T X I 2 H- X for simplicity. We find from 
Q + oc W mag cx TE/T 1 cx £i-2(i-2s)/(7-4s) and 

Q r - ad cx S 2 TV2iJ-l oc S 2+(l- 2s )/(7-4 S ) T l/2 (gee 

Section [3~TT) that the low-/? solutions do not sat- 
isfy the criterion (|57|) . that is, are thermally stable. 
We note that the thermal stability is independent 
of s. 

We also remark the thermal stability of the 
other solutions (ADAF/RIAF, SLE, and LHAF) 
in which the gas pressure is dominant (VFtot ~ 
Wgas cx ST). Equations (HZj) and §T9$ yield 
H cx T 1 / 2 and M cx ST. We find from Q+ cx 
W gas cx ST, Q rad cx tfT^H- 1 cx S 2 , Q ad cx 
MT£, cx ST 2 £ that the ADAF/RIAF solutions 
(Q + ~ Qad) are thermally stable but the SLE 
(Q + ~ Qrad) an d LHAF (Q ad - Q7 ad ) solutions 
are thermally unstable. 

4.3. A Candidate for the Bright/Hard 
State 

The main purpose of this paper is to explain the 
bright/hard state observed during the bright/slow 
transition in the rising phases of the transient out- 
bursts of BHCs. In the low/hard state, the X- 
ray spectrum is described by a hard power law 
with a high energy cutoff at ~ 200 keV. When 
their luminosity exceed ~ O.lTEdd, these systems 



undergo a transition from the low/hard state to 
the bright/hard state. In the bright/hard state, 
the cutoff energy decreases from ~ 200 keV to ~ 
50 keV as the lumino sity increases from ~ O .lTEdd 
to ~ 0.3i E dd (e.g.. iMivakawa et all 120081 ). Be- 
yond the bright /hard state, these systems undergo 
a transition to the high/soft state going through 
the VH/SPL state. 

The ADAF/RIAF solution explains the X-ray 
spectrum in the low/hard state because the elec- 
tron temperature is high (T > 10 9 5 K). How- 
ever, this solution does not exist at the high mass 
accretion rates and does not show the strong 
anti-correlation between the electron tempera- 
ture and the mass accretion rate observed in the 
bright/hard state. The low-/? solution extends 
to such high mass accretion rates. Therefore, 
the disk initially staying in the ADAF/RIAF 
state undergoes transition to the low-/? state 
when the mass accretion rate exceeds M C ,A- On 
the low-/3 branches, the electron temperature is 
low (T c - 10 8 - 10 9 5 K - 10 - 300 keV) and 
strongly anti-correlates with the mass accretion 
rate. These features are consistent with the anti- 
correlation between the luminosity and the en- 
ergy cutoff observed in the bright/hard state. 
Therefore, the low-/3 solution can explain the 
bright/hard state. 

In Section l4~Tj we have discussed the lower limit 
of /? below which the MRI is stabilized therefore 
the low-/3 solutions may not exist. Since /? cx M~ 2 
in the low-/? solutions, this means that the low- 
j3 solution has a maximum mass accretion rate, 
M c ,/9- The mass accretion rate and the electron 
temperature at the lower limit of /? are depicted 
by filled circles in Figure 0] and [5] (see also Table 
[T]). We found that when £ > 0.3 (or equivalently, 
$(t37 = 5r s ) > 4.9$ out ) M C)/ 3 > M C ,A- This in- 
dicates that the disk staying in the ADAF/RIAF 
state undergoes transition to the low-/? disk when 
the mass accretion rate exceeds M Cj a, after that, 
the disk undergoes transition to the optically thick 
disk when the mass accretion rate exceeds M c ,/3- 
In this case, the hard-to-soft transition occurs at 
M c ,/9 (not at M c ,a)- This can correspond to the 
bright / slow transition. 

Here we also remark on the dark/fast tran- 
sition during which the system undergoes tran- 
sition from the low/hard state to the high/soft 
state without going through the bright/hard state 
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and the VH/SPL state. We found that when 
C < 0.3 (or equivalently, $(nj = 5r s ) < 4.9$ ou t) 
M c ,/3 < Mc,A, that is, there is no optically thin, 
thermally stable solution above M C ,A- Therefore, 
the disk staying in the ADAF/RIAF state imme- 
diately undergoes transition to an optically thick 
disk without through the low-/3 disk. This can 
correspond to the dark/fast transition. Accord- 
ingly, we conclude that the bright/slow transition 
occurs when $ has a large value and the dark/fast 
transition occurs when $ has a small value. 

5. Summary 

We have obtained the low-/? solutions for op- 
tically thin, two-temperature accretion disks in- 
corporating the mean azimuthal magnetic fields, 
and concluded that the low-/3 solutions explain the 
bright/hard state observed during the bright/slow 
transition of BHCs. We assumed that the en- 
ergy transfer from ions to electrons occurs through 
Coulomb collisions, and considered bremsstrahhmg 
emission, synchrotron emission, and inverse Comp- 
ton scattering as the radiative cooling processes. 
We prescribed that the Maxwell stress is pro- 
portional to the total (gas and magnetic) pres- 
sure. In order to complete the set of basic 
equations, we specified the radial distribution of 
the magnetic flux advection rate. Accordingly, 
a decrease in temperature can result in an in- 
crease in magnetic pressure under the conserva- 
tion of the magnetic flux in the vertical direc- 
tion. In the low-p solutions, the magnetic heat- 
ing is enhanced by the high magnetic pressure. 
The fraction (1 — <5hcat) of the magnetic heat- 
ing goes into ions and is transferred to electrons 
via Coulomb collisions. The fraction <5hcat of the 
magnetic heating goes into electrons. Eventually, 
the total amount of the magnetic heating goes 
into electrons and balances the radiative cool- 
ing (Compton cooling by bremsstrahlung and/or 
synchrotron photons). The electron temperature 
is lower than that in the ADAF/RIAF solutions 
(T - 10 8 - 10 9 5 K - 10 - 300 keV) and strongly 
anti-correlates with the mass accretion rate in the 
low-/3 solutions. These features are consistent with 
the X-ray spectrum observed in the bright/hard 
state. 



According to lPessah fc Psaltisl (|2005l ). the MRI 
is stabilized for j3 < 2c s /vko- This indicates that 



the low-/3 solutions disappear below this critical 
point. When the magnetic flux advection rate is 
high (C ^ 0.3), the critical mass accretion rate 
of the low-/3 solutions, M c ,/3, is above the maxi- 
mum mass accretion rate of the ADAF/RIAF so- 
lutions, M C ,A- Therefore, the disk initially stay- 
ing the ADAF/RIAF state undergoes transition 
to the low- /8 state when the mass accretion rate 
exceeds M C ,A> after that, the disk undergoes tran- 
sition to the optically thick disk when the mass ac- 
cretion rate exceeds M c .p- This corresponds to the 
bright/slow transition. On the other hand, when 
the magnetic flux advection rate is low (£ < 0.3), 
M c ,/3 is below M c .a- Therefore, the disk initially 
staying the ADAF/RIAF state immediately un- 
dergoes transition to the optically thick disk when 
the mass accretion rate exceeds M C ,A- This can 
correspond to the dark/fast transition. 
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Fig. 1. — Local thermal equilibrium curves of op- 
tically thin, two-temperature accretion disks at 
w = 5r s in the S vs. M/M E dd, 7i(thin line), 
T c (thick line), and (3 plane. The disk parame- 
ters ar e a = 0.05, & = 1, £ e = 0.5, S hcat = 
C = 0.6 (solid), 0.3 (dashed), and 
(dotted). M E dd = L Edd /rj c c 2 = 4nGM/ (rj c K cs c) 
is the Eddington mass accretion rate and rj c = 0.1 
is the energy conversion efficiency The equilib- 
rium curves consist of the ADAF/RIAF, SLE, and 
low-/3 branches. 
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Fig. 2. — Energy balance for ions and elec- 
trons for the same parameters as Figure [1] (solid: 
C = 0.6, dashed: C = 0.3, dotted: C = 0). 
Upper: the ratio of the heat advection to the 
magnetic heating for ions, Q a d,i/(1 - ^hoat)Q + 
(this parameter is referred to as the advec- 
tion factor f. e.g.. iNaravan fc Yil 1 1994 Il995t 



tion tactor J, e.g.. llNlaravan &; Yj liijiMl llVJUat 
lAbramowicz et aLT 19951 Yuan 2001 ). Lower: the 



ratio of the heat advection to the total heating 
(thin), Qad,o/(^hcatQ + + Q IC ), and the fraction of 
the energy transfer via Coulomb collisions to the 
total heating (thick), Q le / (<5heatQ + + Q le ), for elec- 
trons. 
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Fig. 4. — Relations between the electron temper- 
ature and the mass accretion rate for the same 
parameter as Figure [T] Filled circles represent 
critical points at which va — t/CsVkq (or equiv- 
alcntly /3 — 2c s /uko), where va = B v / y/Airpo is 
the toroidal Alfven speed, c s is the sound speed, 
and uk = vjQko is the rotational speed. 



Fig. 3. — Vertically integrated radiative cool- 
ing rates for ( = (bottom), 0.3 (middle), and 
0.6 (top). Solid: bremsstrahhmg (Qbr), dotted: 
bremsstrahlung-Compton (Qbr,c)j dashed: syn- 
chrotron (Q S y), and dash-dotted: synchrotron- 
Compton (Q S y : c)- 
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Fig. 5. — Same as Figure Q] but for a = 0.2. 
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Fig. 6. — Same as Figure [4] but for a = 0.2. 
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Fig. 7. — Dependence of local thermal equilibrium curves on <5h ea t- The disk parameters are a — 0.05, 
£ = 0.5, Co = 0.5, (Wt = (1 + "i^TjTc)- 1 (solid), 0.5 (dashed), 0.2 (dotted), 10~ 3 (dot-dashed), ( = 0.6 
(left) and (right). The exact value of 5hoat is not so important for the low-/? solutions. 
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Fig. 8. — Energy balance for ions and electrons for the same parameters as Figure [7] (left: £ = 0.6, right: 
C = 0). Upper: the ratio of the heat advection to the magnetic heating for ions, Q a d,i/ (1 — £>hcat)<3 + - Lower: 
the ratio of the heat advection to the total heating (thin), Q a d.c/(^hcat<3 + + Q 10 ), and the fraction of the 
energy transfer via Coulomb collisions to the total heating (thick), Q le /(<5hcat<3 + + Q le ), for electrons. 
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Fig. 9. — Local thermal equilibrium curves for the 
same parameters as in Figured] but for £ = —1.0 
(dashed). The equilibrium curves consist of the 
LHAF, SLE, and low-/3 branches. We also plotted 
the thermal equilibrium curves for £j = 1.0 (solid) 
for comparison. 
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Fig. 10. — Local thermal equilibrium curves plotted on the £ - 7] (thin) and T c (thick) plane fc 
£ c = 0.5 (solid), (dashed), —0.5 (dotted), ( — 0.6 (left panel), and (right panel) 
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Table 1 

Values at the lower limit of (3 



a 


c 




M c ^/M Edd 


T c [K] 


0.05 


0.6 


0.27 


1.3 x 10~ 2 


1.7 x 10 9 




0.3 


0.09 


6.4 x 10~ 3 


1.0 x 10 9 







0.03 


4.5 x 10~ 3 


1.3 x 10 s 


0.2 


0.6 


0.26 


1.2 x lO" 1 


1.1 x 10 9 




0.3 


0.08 


4.8 x 10~ 2 


7.4 x 10 s 







0.02 


2.6 x 10^2 


1.3 x 10 s 
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